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The general theory of electromagnetic-fluctuation-induced interactions in dielectric bodies as 
formulated by Dzyaloshinskii, Lifshitz, and Pitaevskii is rewritten as a perturbation theory in terms 
of the spatial contrast in (imaginary) frequency dependent dielectric function. The formulation can 
be used to calculate the Casimir-Lifshitz forces for dielectric objects of arbitrary geometry, as a 
perturbative expansion in the dielectric contrast, and could thus complement the existing theories 
that use perturbation in geometrical features. We find that expansion in dielectric contrast recasts 
the resulting Lifshitz energy into a sum of the different many-body contributions. The limit of 
validity and convergence properties of the perturbation theory is discussed using the example of 
parallel semi-infinite objects for which the exact result is known. 
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I. INTRODUCTION 

Macroscopic material boundaries that interact with 
fluctuating electromagnetic fields experience an induced 
interaction amongst themselves. This was first demon- 
strated for the case of two perfectly conducting paral- 
lel plates by Casimir and was subsequently general- 
ized to take into account frequency dependent dielectric 
properties of the objects by Lifshitz 2\. Relevant exper- 
imental studies were being developed for a long time [3j] 
until recently several high precision measurements of the 
Casimir force were performed ^ . Lifshitz theory and its 
application to various situations such as materials with fi- 
nite conductivity and finite temperature effects has been 
an active area of research in the past few years [11 

The recent trend in miniaturization of mechanical de- 
vices naturally brought up the issue that Casimir-Lifshitz 
interactions need to be taken into consideration when 
small components are at close proximity of each other 
0. For traditional design strategies these forces, which 
could dominate all the others at distances smaller than 
a few hundred nanometers, are to be eliminated. On the 
other hand, one can also imagine using them for novel 
design ideas that could potentially change the way we 
think about designing mechanical systems at that scale 
. For these reasons, it is necessary to develop a better 
understanding of Casimir-Lifshitz interactions when the 
objects involved do not have ideal geometrical shapes. 

This is far from a trivial task, but in recent years 
there have been a number of significant developments 
to this end. These include perturbative approaches for 
geometries that can be considered as slightly deformed as 
compared to some ideal geometries [l^, semiclas- 

sical approaches (Tl| and classical ray optics approxima- 
tions [12 1, multiple scattering and multipole expansions 



[11 [5, [i3|> world-line method [3 and exact 
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numerical diagonalization method [19|, and the recently 
developed numerical Green's function calculation method 
[2^. These different and complementary approaches are 
useful in understanding the subtleties involved in the de- 
pendence of the Casimir-Lifshitz energy on geometrical 
features of the boundaries. 

An interesting result of the Lifshitz theory is that 
the interaction at small separations are effectively de- 
termined by the value of dielectric constants at relatively 
high (imaginary) frequencies where it is not much differ- 
ent from unity [2]. This suggests that a useful comple- 
mentary strategy could be pursued based on expansion 
in dielectric contrast. This approach, which has been the 
subject of a few recent studies [2l|, [23, [H, [H, [2^ , is use- 
ful as it can treat the effect of the geometry of boundaries 
exactly. Here we present a systematic formulation of this 
approach for the calculation of Casimir-Lifshitz energy 
for dielectric objects of arbitrary shape, in the form of 
an expansion in powers of the dielectric contrast in the 
medium (see Fig. [T|). We find that expansion in powers 
of the difference between dielectric constant as compared 
to the background takes on the form of an expansion in 
multi-body contributions to the interaction. We provide 
explicit expressions for each term in the expansion in the 
form of convolutions of a tensorial kernel. We show that 
a resummation of the expansion can help augment the 
convergence properties of the series, and make it appli- 
cable to a wider range of dielectric properties. 

The rest of the paper is organized as follows. Section 
lUlays out the general formalism used for calculating the 
Casimir-Lifshitz interaction, which is applicable in gen- 
eral to metals and dielectrics. In Sec. IIIIl the dielectric 
contrast perturbation theory is developed based on the 
formalism of Sec. |TT] and an explicit expression is ob- 
tained for each term in the expansion. In Sec. IIVI the 
series obtained in Sec. Illll is resummed using a decompo- 
sition of the kernel involved in the general expression for 
the energy. Sect ion fVl gives the derivation of the Casimir- 
Polder energy using the formalism developed in Sec. llVi 
as an example. The convergence properties of the series 
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FIG. 1: Schematics of the dielectric function profile, and the 
many-body contributions to the electromagnetic-fluctuation 
forces. 



obtained as well as the nature of the divergent contri- 
butions in the theory are discussed in Sec. IVII Section 
I VIII is devoted to a specific class of geometries where two 
nearly parallel semi-infinite dielectric objects are placed 
in front of each other. Finally, Sec. IVIIII concludes the 
paper with some discussions. 



we choose to describe in the temporal gauge, where the 
electrostatic potential vanishes identically, namely = 0. 

We start from the thermal (finite temperature) Green's 
function of the electromagnetic field (Cs; r, r') in imag- 
inary frequency, where Cs — 2T:sk^T /h is the Matsubara 
frequency with ^rT being the thermal energy and s a 
positive integer 26]. In the medium described above, the 
thermal Green function is the solution to the following 
equation 



By introducing the operator 



I?jfe(G;r,r') = <5,fc,53(r-r'). 

(1) 



/Cy(Cs;r,r') 
we can write Eq. ([T]) as 



(2) 



which means that as an operator we have 



/C.,(C.;r,r') = I?,^^(Cs;r,r'). 



(4) 



II. GENERAL FORMALISM 

Let us assume that we have an assortment of dielec- 
tric objects in space with arbitrary shapes and frequency 
dependent dielectric properties, as sketched in Fig. [T] 
This medium can be described using a frequency- and 
space-dependent dielectric function e(w,r). We consider 
fiuctuating electromagnetic fields in this medium, which 



Now imagine a process in which we start from the empty 
space and introduce the dielectric objects into the space 
in a perturbative manner, similar to the charging pro- 
cess of a capacitor or constructing a charge distribution 
by bringing infinitesimal charge elements from infinity to 
assemble the distribution. Using standard diagrammatic 
techniques f27| , we can show that the introduction of the 
dielectric objects causes the Helmholtz free energy of the 
system to change according to the following formula (26l | 



ST = 5To + kBTY, /d'rid3r2l?y(G;ri,r2)(57',,(Cs;r2,ri), (5) 



s=0 

where 

a2 



SV,,iO;r,r') = ^^..^^(r - r')5e(*Cs, r), (6) 

is the Polarization operator, and the primed summation means that the s = term has an extra factor of i. In the 
diagrammatics, the Polarization operator is formally related to the Green function via the Dyson equation 



where pj^-* ((s; r, r') denotes the Green function before of new material. Equation ([7]) can be solved to yield 
the change in the dielectric profile due to the introduction 

P„.i(C.;r,r') =P,(°^"'(C.;r,r') + m,(C;r,r'). (8) 
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Equations ([2]) and ^ lead us to the following observation 
<57'.,(C.;r,r')=<5/C,,(G;r,r'). (9) 
which helps us to write Eq. ([5]) in the form of 



SJ' = 5Tn + fceT^ tr [JC-^5K] 



(10) 



Equation (|T0|) can be formally integrated to yield the fol- 
lowing expression for the contribution to the Helmholtz 
free energy due to Casimir-Lifshitz (CL) interactions: 



^CL = fcBr^trln[/C,,(Cs;r,r')] 



(11) 



s=0 



At (effectively) low temperatures, we can convert the 
summation over s into an integration by changing As = 1 
into ds — hdC, / {2nk'QT) . In this case, we find the 
Casimir-Lifshitz energy as 

i?CL = ftj|"^trln[^,,(C;r,r')]. (12) 

This could be a relatively simple starting point for the 
calculation of the Casimir-Lifshitz energy in any system. 



III. DIELECTRIC CONTRAST 
PERTURBATION THEORY 

We can now construct a systematic perturbation the- 
ory scheme based on the above definition of the Casimir- 
Lifshitz energy. The starting point is to write the dielec- 
tric function profile as 



e«,r) = l + 5e(iC,r), 



(13) 



which can be used to decompose the kernel /Cy into a 
part that corresponds to the empty space and a pertur- 
bation that entails the dielectric inhomogeneity profile. 
In Fourier space, this reads 



(C; q, q') = (C, q) {'2irf5\^ + q') + ((; q, q'), 
where 



and 



^0,jj(C,q) = -^^t] + (fS^j - q^qJ, 



5/Cjj(C;q,q') = ■^(5ij(5e(iC,q + q') 



(14) 
(15) 



(16) 



This decomposition can now be used to construct the 
perturbation theory. 

The expressions for the Casimir-Lifshitz energy [Eqs. 
(fTTj) and mi)] involve trln [/C], which can be written as 
a perturbative series by using 



trln[A:] = trln[A:o] +trln[X + /Coi(5/C] 



trln[/Co] + ^ 



tr 



(/C^M/C)' 



where I is the identity tensor and 

^o;l,(C,q)- 











c2 





(17) 



(18) 



Putting in the explicit forms for /Cg ^ and (5/C, we can 
write the explicit form for the trace as 



tr[(/Coi(5/C)"i 



(27r)3 (27r)3 



(19) 

The above expression contains the geometric information about the arrangement of the dielectric objects through the 
Fourier transform of the dielectric function profile. We can now rewrite Eq. (|19p in real space, and find the following 
series expression for the Casimir-Lifshitz energy of any heterogeneous dielectric medium 



dC ^ (-1) 



27r 



E 



ri-l 



d^n . . .d3r„ g^^^^ir, - r2) • • • ^t,^(r„ - n) 5e«,ri) • • • 5e(zC,r„), (20) 



where the Green function defined as &j^,(r) = ^ICqIA^jT) reads 



2 p-C/c 



0,ij \ 

c c 

1 + 



nr. 

r 



2 



c c 
1 + 3— -f 3 



+ x'5,,^'(r). 



(21) 



in real space. The tensorial kernel Q^, , has the struc- ture of the electric field of a radiating dipole in imagi- 
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nary frequency, including the delta function contribution 
that ensures appropriate behavior of the expression in 
the near-field [23| . This kernel has been introduced some 
time ago in connection with van der Waals interactions 

The main result of Eq. (j20[) [and a re-summed version 
of it given in Eq. below] has a number of interest- 
ing characteristics. First, We find that an expansion in 
powers of Se automatically turns into a summation of in- 
tegrated contributions of n-body interactions. Moreover, 
all of the n-body interaction terms have simple explicit 
expressions in terms of a single fundamental kernel that 
mediates the two-body part of the interaction. Finally, 
we find a closed form expression for the Lifshitz energy 
for any geometrical arrangement of dielectric bodies in 
terms of quadratures, which is amenable to simple di- 
agrammatic rules and can hence be easily adapted for 
numerical computations at any given order of the per- 
turbation theory. 



IV. CLAUSIUS-MOSSOTTI RESUMMATION OF 
THE PERTURBATION THEORY 

The perturbation theory developed above uses the spa- 
tial contrast in the dielectric function as expansion pa- 
rameter and this in general may not be a suitably con- 
trolled expansion parameter, especially at zero frequency 
where the dielectric contrast could be considerably larger 
than unity (see Sec. IVII for more discussion). The for- 
mulation can be augmented by resumming the pertur- 
bation theory such that it is organized as an expansion 
in powers of the combination i^^g^^^^ = instead 
of 6e. The new expansion parameter, which reminds us 
of the Clausius-Mossotti equation for molecular polariz- 
ability ^28], is systematically smaller than the dielectric 
contrast itself, and is finite even for real metals at zero 
frequency (where the dielectric contrast diverges). 

The key to this remedy lies in the Dirac delta function 
term in the expression of the kernel in Eq. (|2ip . This 
suggests a decomposition of the form 



/Co 



(22) 



to be used in the formulation, where the operator A is 
defined as 



Aj(C,q) 

in Fourier space, and as 

^2 g-Cr/c 



>2 „ 

2^i5y + ^QiQ] - q Si. 



(23) 



Aj(C,r) 



1 + ^ 



1 + 3^ + 37^ 



2 



(24) 



in position space. Putting Eq. p2|) in Eq. (|17p . we find 



tr In [K] ^ tr In [/Co] + tr In [X + K^^SIC] , 
tr In [/Co] + tr In {l+ —XSK. 



tr In [/Co] + tr In 



3(2 
^2 



rA5K 



r(5/C 



-trln 



Defining a new operator 
SB=-^SIC ( 



3C 



:6IC 



(25) 



(26) 



which in position space has the following explicit form 



l + i,5e(zC,r) 



<5%(C;r,r') = 
we can rewrite Eq. ([25]) as 
trln[/C] = trln [/Co] +trln 
= tr In [/Co] + tr In 

Z — / ri 



S,,SHr-v'), (27) 



3(2 



r(5/C 



trln[I + ASB] 



3(2 
tr [{ASBy 



SIC 



(28) 



This form can now be used to recast the perturbative 
expansion for the Casimir-Lifshitz energy as a systematic 
expansion in powers of the new perturbation parameter, 
which is always less than unity. We thus find 



CL 



= h 



dC y> (-1)"'^ 



27r 



(TYi ■ ■ ■ (TVn 



X Aii2(i"i - r2) • ■ •-4i„,;i(r„ - ri) 



(Se(iC,ri) 



l + \Se{iC,r{) 



'5e«,r„) 



1-K i(5e«,r„) 



,(29) 



in real space. In this result, we have neglected two re- 
maining terms in Eq. (|28[) . The first one is the trivial 
term trln [/Cq] that corresponds to the self energy of vac- 
uum and could be easily eliminated as it does not depend 
on any of the physical parameters involved. The second 
term in Eq. (j28p is a singular contribution, which has 
the explicit form of 



d^r In 



(2^ 



(30) 



This term does depend on the geometry of the dielectric 
objects, and need to be subtracted off for any given ge- 
ometry and dielectric configuration. We will comment on 
the general issues related to such divergent contributions 
as well as the convergence properties of the series in the 
following Section. 
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V. EXAMPLE: CASIMIR-POLDER 
INTERACTION 

To see how the formalism works, let us use it to cal- 
culate the Casimir-Polder (CP) interaction between two 
objects at long separations. Consider a sphere of volume 
vi and dielectric function ei(iC) located at the origin, 
and a sphere of volume V2 and dielectric function e2(*C) 
located at a separation R from the first one. To calcu- 
late the Casimir-Polder energy, we can use Eq. ((29|) in 
conjunction with 



Ignoring the frequency dependence of the polarizabilities, 
we find 



«2(53(r-R). (31) 



This is a singular form for the profile and is only valid 
in the limit that the separation is much larger than the 
typical size of the objects. Therefore, to be consistent 
with this singular limit we should only keep the leading 
order term in the size of the objects. A typical term 
in Eq. ([^5]) with this dielectric profile, which does not 
contribute to the self energies of the two objects, would 
(symbolically) look like 



l + |^ei«). 



+ hSe2iiO 



which involves multiple powers of the quantity 

/ = «A,(0). (33) 

Using the proper definition of / within our singular de- 
scription of the objects, we find 



/ 



lim 



An 



A,(r) 



0. 



(34) 



Therefore, to the leading order, we find 



27r 



1 + 4-561(^0 



V2Se2iiC) 
l + kSe2{iC) 



X Aj(R)Az(-R). (35) 
Using the following definition for dynamic polarizability 



^i'^ei«) 



(36) 



and Eq. 



Ecp — 



we can rewrite Eq. (|35p as 

h 



3 + 6 — 
c 



dC ai{i()a2iiO e 

C'i?2 (3^3 ^4^4 



.(37) 



23 he 
Ecp = ~ — -=^OL\a2, 

4:71 R' 



(38) 



which is the celebrated result for Casimir-Polder interac- 
tion ^3d\ . The above results could also be obtained using 
Eq. ([20| in which case the corresponding /-factors would 
create nonvanishing constant values that would need to 
be added up to generate the Clausius-Mossotti form of 
the polarizabilities. 

While the limiting form at longest separations can be 
obtained by this simple treatment, the correction terms 
(that are appreciable at closer separations) will systemat- 
ically be produced from to the /-factors that technically 
speaking correspond to higher multipole contributions to 
the Casimir-Polder energy, as well as many body contri- 
butions. To keep track of these corrections in a system- 
atic way requires a multipole expansion type approach of 
the type developed in Refs. [l5i.Tl^. IT7| . Finally, we note 
that the simple treatment given above corresponds to 
spheres where the symmetry of the object easily guaran- 
tees that Eq. ([M]) holds true. For more complicated ge- 
ometries, the question of orientation as well as the shape 
of the object complicates matters more [l^ [l3| and a 
more systematic approach is needed. 



^+1 (— R-)A^+iz„+2 (O2) • • • A„ii (R), 

(32) 



VI. CONVERGENCE AND REGULARIZATION 
OF THE PERTURBATION THEORY 

A. When is the series expansion convergent? 

Because the formulation is constructed as a perturba- 
tion theory, which can perhaps complement other pertur- 
bation theories for Casiniir-Lifshitz interaction between 
objects of arbitrary geometry, we need to address the 
convergence properties of the expansion. To examine the 
convergence property of Eq. (|20p . we can use a specific 
example for which the exact Casimir-Lifshitz energy is 
known, and compare it with the perturbative approxi- 
mations at given orders. We consider the case of two 
identical semi-infinite dielectric objects with flat bound- 
aries that are placed parallel to each other at a separation 
H . The exact expression for the energy per unit area for 
this problem is known to be 



^Lif — 



In 



In 



dC c 



{s + pf 



dp p 



-2pCH/c 



-2pCH/c 



(39) 



where s = y^e — 1 + p'^. We assume a simple form of 



s«) = 1 + 



(40) 
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FIG. 2: Comparison of the expansion in Eq. (|20p up to the 
2nd, 4th, 6th, and 8th order, with the exact Lifshitz result 
for a pair of parallel semi-infinite dielectrics, as a function of 
their separation H in units of the plasma wavelength Xp. The 
dielectric constant of both objects is chosen as Eq. (|40|) with 
(a) LJo = 1.4 LJp, (b) ujo = 1.0 Wp, and (c) ujo ~ 0.8 uip. 



FIG. 3: Comparison of the expansion in Eq. (I29p (the so- 
called Clausius-Mossotti resummation) up to the 2nd, 4th, 
6th, and 8th order, with the exact Lifshitz result for a pair 
of parallel semi-infinite dielectrics, as a function of their sep- 
aration H in units of the plasma wavelength Xp. The dielec- 
tric constant of both objects is chosen as Eq. (|40p with (a) 
ujo = 1.4 u!p, (b) UJO = 0.8 tjp, and (c) ojq = 0.4 uip. 



for the dielectric constant (in imaginary frequency), 
where ujp represents the plasma frequency, from which 
the plasma wavelength Xp = 27rc/a;p can be extracted. 
We have checked that a dissipative term of the form 7^ 
in the denominator of Eq. (j40p would not alter are re- 
sults significantly for realistic values of 7, therefore we 
have neglected this term for simplicity of our presenta- 
tion. 

In Fig. [21 the ratio between the energy as calculated 
from Eq. (PO)) and the exact Lifshitz result is shown as 
a function of the separation, up to the second, fourth. 



sixth, and eighth order in perturbation theory. The re- 
sults clearly show a crossover between two asymptotic 
regimes near H = Ap/(27r) [3l|. Figure [2^ corresponds 
to Wo — 1.4 LOp and shows a rapid convergence. This can 
be understood from the fact that even at zero frequency 
where the dielectric contrast has its largest value, this 
example yields (5e(0) — 0.5 that is smaller than unity. 
Figure [^b shows that when = 1-0 the series is still 
convergent (despite 5e(0) = 1), although not as rapidly 
as in the example of Fig. [2^. For loq = 0.8 Wp, the 
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series is divergent, as Fig. [2]: shows. Therefore, the se- 
ries in Eq. (|20p appears to be rapidly convergent when 
Se{0) < 1, and the convergence is considerably more ef- 
ficient for H < Ap/(27r) [as compared to H > Ap/(27r)], 
especially for ^e(O) < 1. 

The resummation of the perturbation theory discussed 
in Sec. IIVI improves the convergence property of the se- 
ries as it effectively replaces Se (that can have any value) 
by — ' "^hich is bound to be between and 

3 even for real metals at zero frequency. More explicitly, 
using Eq. (j40p one can see that the expansion parameter 
at zero frequency, which is the parameter that controls 
the convergence of the series at large separations, changes 
from 6e{0) = Wp/wp ^^^/(Swq + LUp). This suggests 
that at large separations, the series [in Eq. ([H])] should 
now be convergent for too > 0.8 ujp (so that the expan- 
sion parameter is smaller than unity). Figure [3] shows the 
convergence property of a few examples using Eq. (j29p 
instead of Eq. (|20p . For luq = 1.4 tUp, a much more rapid 
convergence is observed as shown in Fig. [3^, while the 
previously divergent case of = 0.8 ujp now shows a 
good convergence, as seen in Fig. [Sb and in agreement 
with the argument above. As Fig. [St shows, even for 
ujo = 0.4 ujp one still observes a good convergence despite 
the fact that the zero frequency expansion parameter is 
equal to 2. 

Figure [4^ shows the result of the expansion for uiq = 
0.3 LOp. Interestingly, it appears that in this case the 
series is (at least up the eighth order in perturbation) 
convergent for H < Ap/(27r), while it clearly diverges for 
H > Ap/(27r). One can roughly say that the appropriate 
parameter that controls the convergence at small sepa- 
rations is the expansion parameter at — luq, which is 
3ujp/{6ujQ + ujp). For ujo — 0.3 ujp, this parameter is ~ 2, 
which is consistent with the value of the zero frequency 
expansion parameter in the borderline convergence case 
at large separations obtained for wq = 0.4 ujp. Finally, for 
Wo = 0.1 Wp, the series is divergent, as can be seen in Fig. 

It is interesting that the criterion for convergence 
seems to be that the relevant expansion parameter — i.e. 
the parameter at the relevant frequency — needs to be 
smaller than 2. 



B. Divergencies and Regularization 

A general issue with the calculation of the Casimir- 
Lifshitz energy is the appearance of divergent contribu- 
tions, like in any quantum field theory calculation. In 
the case of a system with arbitrary geometry, it is not 
clear how these divergencies depend on the details of the 
geometry, so that they could be identified and dealt with 
in a systematic way. The present formulation also suffers 
from the presence of divergent terms, but similar to any 
perturbative field theory with a manifest recipe for con- 
structing each term, these divergent contributions could 
be systematically identified and regularized order by or- 
der in perturbation. 
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FIG. 4: Comparison of the expansion in Eq. (I29p (the so- 
called Clausius-Mossotti resummation) up to the 2nd, 4th, 
6th, and 8th order, with the exact Lifshitz result for a pair 
of parallel semi-infinite dielectrics, as a function of their sep- 
aration H in units of the plasma wavelength Xp. The dielec- 
tric constant of both objects is chosen as Eq. (|40p with (a) 
LJo = 0.3 Ljp (corresponding to silicon) and (b) uio = 0.1 uip. 



As pointed out by Barton [21|, it is important to ex- 
amine the nature of the divergent contributions, and for 
example determine whether they are controlled by the 
minimum possible distance between atoms and molecules 
(i.e. cutoff on wavevector), transparency of the materials 
at high frequency (i.e. plasma frequency as cutoff) or ul- 
traviolet frequency cutoff in vacuum. This point has not 
always been carefully dealt with in the literature of the 
Casimir-Lifshitz interactions, partly because the largely 
used perfect conductor limit already blurs this distinc- 
tion at the outset when it assumes the plasma frequency 
is infinite. 

In our formulation, the divergencies originating from 
lack of molecular excluded volume in the theory will be 
regularized using a high wavevector cutoff in the q inte- 
grals in Eq. ([TO]) or a short distance cutoff in the position 
integrals in Eqs. ([20]) and ((29l) . To see how this works 
let us go to Eq. (|19p and look at the q integrations. We 
have 



tr[(/Co-M/C)"] = 
X 5e{iC 



(1) 



rf3q(i) 
(27r)3 
fq(2)) 



d3q(«) 
"(2^ 



•••,q("^) 

[(!)). (41) 



where 



M(q' 



(1) 



q(i).q(2))...(q(n).q(i)y / 

(42) 



Jn)2 



Since the (Fourier-space) dielectric contrast profiles only 
depend on the differences between the wavevectors, we 
can change the integration variables so that one inde- 
pendent wavevector can be integrated out. Using = 



Q where Q 



tr[(/Co^i5/C) 



.(1) 



(27r)3 
X M(qi + Q, 

X Se{iC, -qi + q2) 

where q„ = — qi — • • • — q„_i. The integral 



q(")), we find 
(27r)3 (27r)3 

, q™ + Q) 

<5e(zC,-q„ + qi)j(43) 



d^Q 

(27r)3 



M(qi + Q,- 



Q), 



can now be performed independently, and looking at 
the definition of the kernel M in Eq. ([^^ we can see 
that it diverges at high values of Q where M ^ 1. 
This divergent contribution can be extracted by writing 
M = Af — 1 + 1 and separating the 1 factor, which leads to 
a regularized contribution proportional to / — 
The remaining q integrations for this divergent contribu- 
tion can be carried out, and the series can be summed 
up to yield an overall singular contribution of the form 



dC 

2 



^ J d^r ln[e«,r)] 



(44) 



which corresponds to volume terms and can be system- 
atically isolated for any geometry and dielectric configu- 
ration. This result entails a similar volume contribution 
calculated in Ref. [2l| at the second order in Se. Other 
divergent contributions in the present formulation can be 
dealt with in a similar manner to the above example. 

Finally, we note that using a very similar calculation 
to the one presented in Sec. |V]for the Casimir-Polder in- 
teraction, one can show that for an arbitrary assortment 
as depicted in Fig. [1] the Casimir-Lifshitz energy calcu- 
lated from the present formalism in the limit where the 
objects are far from each other only contains divergent 
contributions in the self energies of these objects and all 
of the terms that depend on the distances — including the 
many-body terms — are finite. 



VII. PARALLEL SEMI-INFINITE OBJECTS 

Due to its importance, we now focus our attention on 
the specific arrangement shown in Fig. [51 where two 




FIG. 5: Scliematics of two semi-infinite dielectric bodies with 
boundaries of non-ideal shape. 



nearly parallel semi-infinite dielectric bodies with irregu- 
larly shaped boundaries are placed next to each other at 
a mean separation H. We can write down the dielectric 
function profile in space as 



e2(iC), H + h2ix) < z < +00, 
e«,r) = <(l, /ii(x) < z < i/ + /i2(x), 

ei(iC), — oo < z < hi{x), 
which in Fourier space reads 

5f«,q) - — / d^xe'i-L '' 
Qz J 



(45) 



S€2 e 



iq. [H+h2{x)] 



Sei e 



iq^hi (x) 



(46) 



When the parameters are such that the leading contribu- 
tion in Eq. ((29|) comes from the second order term (see 
Sec. IVip , we can simplify the expression for the Casimir- 
Lifshitz interaction and write it in a closed form that can 
be readily used in studies of various geometrical effects. 
A similar strategy has been discussed in Ref. [2l[ , where 
a large variety of other geometries has been considered. 
Putting in the dielectric function profile of Eq. P5)) . we 
find 



^2 



1287r3c4 7o 



dC C 



ei(*C)-l 



e2(*C)-l 



e2(*C) 



d^xd^x' w{-[x-x%-[H + /i2(x) - /ii(x')] 
c c 



(47) 



where 



W(y, h)^8T (O, 2 + ^ 



s3/2 

4 (1 - h^s) 



[y^^K^sY {y^ + h^s)^^ {y'^ + h^s) 



(48) 



and r{a,z) = dt t""'^ e~* is the incomplete gamma 
function. Note that at this order, the Lifshitz energy is 
pairwise additive, and that the effect of geometry has 
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been taken into account exactly for any arbitrary profile. 
In the above result, we have kept the frequency depen- 
dence of the dielectric functions as well as the geometry 
of the boundaries arbitrary for generality of the presen- 
tation. The expression in Eq. (|47)) can be considerably 
simplified for ft.i(x) = 0: 



(49) 



in terms of the original Lifshitz result for the energy per 
unit area of flat boundaries (within the same scheme of 
Clausius-Mossotti approximation) 0] 



S2{H) = 



647r2c2 



ei«) + 2 



e2(*C) - 1 



(50) 



with 



-2u 



Ciu) = 'i{u-l)Ei(2u) + ^r-n+2u + u^ -iu-"), (51) 

and the exponential integral function defined as £'„(z) = 
dt e~^*/^"- This simplification is a general feature 
that is present for any pairwise additive interactions, as 
has been shown in Ref. @. 



VIII. DISCUSSION 

When dealing with complicated quantum field theories, 
it is always helpful to try and formulate complementary 
perturbative schemes so that the specific point of inter- 
est, which is usually out of reach, is approached from 
different directions. This could potentially provide com- 
plementary information that could be compiled to yield 
an improved insight into the properties of the theory. A 
good example of this synergy is the perturbation theory 
and 1/n expansion of the 0{n) model of the (j)'^ quantum 
field theory [s^. A similar strategy has been the under- 
lying motivation for the work presented here: expansion 
in dielectric contrast allows us to formulate a perturba- 
tive scheme for calculating the Casimir-Lifshitz interac- 
tion between object with arbitrary geometry keeping the 
effect of the geometry exact. The motivation behind this 



alternative approach comes from the fact that at dis- 
tances smaller than the plasma wavelength the effective 
dielectric contrast that determines the Casimir-Lifshitz 
interaction corresponds to the high frequency limit of 
Eq. (|40p . which is effectively smaller than the equiva- 
lent value for the large distance regime. This view is 
confirmed by the plots in Figs. [H [31 and 21 which show 
that the series expansion in dielectric contrast converges 
more rapidly at short distances. Interestingly, Fig. [4^ 
suggests that it is possible to have a convergence at dis- 
tances smaller than the plasma wavelength — typically a 
few hundred nanometers — while the same series diverges 
at larger distances. This marginal case is observed for 
0Jn = 0.3 ujp, which coincidentally corresponds to silicon 
[33| that is a common choice for fabrication of small me- 
chanical components with fine geometrical features. 

The perturbation theory seems to have good conver- 
gence property for materials that have a dielectric con- 
trast of ~ 2, or a dielectric constant of ~ 3, at the rele- 
vant frequency that is zero for large distance asymptotics 
and ujQ for short distance regime. While this is somewhat 
restrictive, as for example it does not include real met- 
als, it is not too far from being able to include some 
commonly used dielectric materials, such as silicon for 
example. Note that this complementary approach keeps 
the effect of the geometrical features exact, which is a 
significant improvement compared with complementary 
theories that apply to objects with small deformations. 
Finally, we note that it might be possible to use Borel 
summation method ^32,] (or other equivalent techniques) 
to improve the convergence of the series expansion in the 
strong coupling limit of dielectric contrast. 

In conclusion, we have developed a perturbative 
scheme for calculating the Casimir-Lifshitz interaction 
between objects of arbitrary geometry. Explicit expres- 
sions are provided for each term in the perturbation the- 
ory, as multiple integrals over the bodies of the inter- 
acting objects. This method could in principle be used 
to calculate the interactions at any order using standard 
numerical integration techniques. 
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